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Abstract 



A new Stokes inversion technique that greatly reduces the effect of the spatial point spread function of the telescope is used to constrain 
the physical properties of umbral dots (UDs). The depth-dependent inversion of the Stokes parameters from a sunspot umbra recorded 
with Hinode SOT/SP revealed significant temperature enhancements and magnetic field weakenings in the core of the UDs in deep 
photospheric layers. Additionally, we found upflows of around 960 m/s in peripheral UDs (i.e., UDs close to the penumbra) and 
« 600 m/s in central UDs. For the first time, we also detected systematic downflows for distances larger than 200 km from the UD 
center that balance the upflowing mass flux. In the upper photosphere, we found almost no difference between the UDs and their 
diffuse umbral background. 
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1. Introduction 

Umbral dots (UDs) are small brightness enhanc ements in 
sunspo t umbrae or pores and were first detected by IChevalierl 
d!916l) . The strong vertical magneti c field in umbrae suppresses 
the energy transport by convection (iBiermannl Il94ll) . but some 
form of remaining heat t ransport is needed to explain the ob - 
served umbral brightness (lAdiabshirzadeh & KoutchmvT,ll983l) . 
Magnetoconvection in umbral fine structure, such as UDs and 
light bridges, is thought to be the main contributor to th e energy 
transp o rt in the umbra (IWeissL I2002), see reviews by ISolankil 
d2003h . lSobotkal(l2006l) . and lBorrero & Ichimotol(l201 ll) . 

Progress in the physical understanding of umbral dots was 
made wi th numerical simulat i ons of 3D radiative magnetocon- 
vection (ISchussler & Voglerl l2006t iBharti et all l2010t) . Most 
of the simulated UDs have a horizontally elongated shape and 
show a central dark lane in their bolometric intensity images. In 
the deepest photospheric layers, the inner parts of UDs exhibit 
magnetic-field weakenings and upflow velocities. The simulated 
UDs are surrounded by downflows that are often concentrated 
in narrow downflow channe ls at the endpoints of the dark lanes 
( ISchussler & Voglerl |200oT ). Higher up in the photosphere, the 
UDs in the simulations do not differ significantly from the dif- 
fuse background. 

Considerable efforts on the observational side were made to 
test these theoretical predic tions. Dark lanes in side UDs were 
found in the observations of Bharti et aT| (l2007b with the 50-cm 
Hinode telescope and by [Rimmeld d2008l) , who observed with 
the 76-cm Dunn Solar Telescope. However. lLouis et al ] d2012h 
analyzed straylight-corrected Hinode/BFI data and did not find 
dark lanes in their observed UDs, which leaves room for doubt 
whether the observed phenomena are really identical with the 
synthetic ones. The UDs described in [B harti et al. I d2007h dif- 
fer from those reported in lSchiissler & Voglerl ( 20061 in that the 



area of the observed features is an order of magnitude larger; 
possibly they are the remains of a decayed light bridge. 

More important than the dark lanes are th e flows, since they 
are cent ral to the convective nature of the UDs. lRiethmuller et al.l 
(2008a), using inversions of Hinode/SP data, discovered upflows 
in the deep layers of peripheral UDs (PUDs) but not in central 
UDs (CUDs), wh ile downflows were not detected. Subsequently, 
lOrtiz et al.l (1201 Oh studied a small pore recorded with the CRISP 
instrument of the 1-m Swedish Solar Telescope and found irreg- 
ular and diffuse downflows in the range 500- 1 00 m/s for a small 
set of five UDs. In contrast, in their recent study, IWatanabe et al.l 
d2012h analyzed a larger set of 339 UDs, also observed with 
CRISP, and found significant UD upflows, but no systematic 
downflow signals. Thus, the existence of downflows in or around 
UDs remains uncertain, so that the fate of the material flowing up 
in UDs is unclear. The depth-dependen t invers ions of full Stokes 
profiles deriv ed inlSocas-Navarro et aLl(l2004b and later at higher 
resolution in iRiethmiiller et al.1 (l2008al) revealed a temperature 
enhancement and a field weakening for the UDs compared to 
the nearby umbral background, which both were strongest in the 
deepest observed layers. 

Since the observational picture is inhomogeneous, there is a 
need for a more detailed UD study for which high spatial and 
spectral resolution is of utmost imp ortance. In t his wo rk, the im- 
proved Stokes inversi on method of Ivan Noortl (1201 2l) is applied 
to Hinode/SP data (see lvan Noort et al.Ll2013l) . This so-called 2D 
inversion method allows the depth-dependent structure to be ob- 
tained basically as it would be in the absence of the telescope's 
point spread function (PSF). 



2. Observation, data reduction, and analysis 

The data we analyzed in this study were recorded from 12:43 
to 1 3:00 UT on 2007 January 5 with the spectropolarimeter 
(SP, iLites et all l200ll) of the Solar Optical Telescope (SOT, 
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Figure 1. Stokes / continuum intensity of the Hinode/SP map of 
a sunspot umbra of NOAA AR 10933. The original data are plot- 
ted in the top panel. The Stokes / continuum resulting from the 
2D inversion is shown in the bottom panel. The intensity is nor- 
malized to the mean quiet-Sun intensity Tqs. The outer contour 
line in the bottom panel indicates the edge of the umbra as re- 
trieved from the magnetic field inclination map (see main text), 
the inner contour line separates central from peripheral umbral 
dots (UDs). Four typical UDs are marked by circles and letters. 

Tsuneta et all 120081) on the Hinode spacecraft dKosugi et all 
2001 " The SP was operated in its normal map mode, i.e., the 
integration time per slit position was 4.8 s, resulting in a noise 
level of 10~ 3 (in units of the continuum intensity). The sam- 
pling along the slit, the slit width, and the scanning step size 
were 0"16, the spectral sampling in the considered range from 
6300.89 to 6303.26 A was 21 mA pixel" 1 . The center of the ob- 
served umbra was located very close to the disk center, at a he- 
liocentric angle of 2.6°. The full Stokes profiles were corrected 
for dark current as well as flat-field effects and calibrated with 
the sp_prep routine of the SolarSoft package. Part of the cali- 
brated Stokes / continuum intensity map obtained with Hinode 
SP is shown in the upper panel of Fig. Q] The original FOV is 
much larger and contains quiet-Sun regions that are used for the 
intensity normalization. 

Under the assumption of local thermodynamic equilibrium, 
the Stokes profiles of the Fe i 6301.5 A and 6302.5 A lines were 
inverted by applying the versio n of the SPINOR inve rsion code 
dFrutigeri 120(3(31 iFrutiger et all l2000b extended by Ivan Noortl 
(2012). In this version of the code, the instrumental effects re- 



sponsible for the spectral and spatial degradation of the observa- 
tional data are taken into account, so that the inverted parame- 
ters correspond to spatially deconvolved values (but without the 
added noise that deconvolution generally introduces). The obser- 
vational data are spatially upsampled by a factor of two, so that 
the input and output data of the SPINO R inversion have a sam - 
pling of 0''0 8 per pixel ( for de tails, see Ivan Noort et al 1 12013b . 
According to Ivan Noortl (|20 1 2b . the spatial PSF used for the in- 
version considers the 0.5 m clear aperture of the SOT, the central 
obscu ration, the spider dTsuneta et all 120081 : ISuematsu et all 
|2008[), and a defocus of 0. 1 waves, even though the focus posi- 
tion of our data set was not accurately known. Height-dependent 
temperature, LOS velocity, magnetic field strength, field inclina- 
tion, field azimuth, and micro-turbulence are determined at three 
log T500 nodes: -2.5, -0 .9, and 0. More details o f th e inversion of 
this sp ot are provided in Ivan Noort et al 1 (l2013h and lTiwari et all 
(f2Q13h . 

A continuum map obtained from the best-fit Stokes / pro- 
files of the 2D inversion result can be seen in the bottom panel 
of Fig. Q] Since the deconvolution of the data with the theoret- 
ical spatial PSF is now indirectly part of the inversion process, 
the contrast is significantly enhanced and UDs can be identi- 
fied much more clearly than in the original data. Hence, the 
continuum map in the bottom panel of Fig. [TJ was used for a 
manual detection of the location of the most prominent 67 UDs. 
Misidentification of brightness features that are separated from 
the penumbra in the continuum image but are still connected 
to the penumbra in the magnetic field inclination map were ex- 
cluded since we defined the boundary of the umbra by thresh- 
olding the lowpass-filtered inclination map (7x7 pixels) at 40 
degrees and corrected the boundary found in this way manually 
in a few doubtful cases where the penetration of long and narrow 
penumbral filaments led to a wrong result. The umbral boundary 
is shown as the outer contour line in the bottom panel of Fig. Q] 
Furthermore, we divided the set of 67 UDs into 23 central UDs 
(CUDs) and 44 peripheral UDs (PUDs). The inner contour line 
in the bottom panel of Fig.[T]separates the PUDs from the CUDs. 
The criterion used is simply the distance to the umbra-penumbra 
boundary. CUDs are > 1000 km away from this boundary, PUDs 
< 1000 km. 

Once the locations of the UDs' centers were known, the UD 
boundaries were determined from the continuum map by a mul - 
tilevel tracking (MLT) algorithm (see iBovelet & Wiehii 1200 ll) . 
using 25 equidistant intensity levels. The resulting contiguous 
MLT structures were then cut at 50% of the local min-max in- 
tensity range, which was taken as the UD boundary. A detailed 
descript ion of the use of the MLT algorithm for isolating UDs is 
given in lRiefhmuller et all d2008bl) . 

The knowledge of the UD boundaries allowed us to av- 
erage UD properties over all pixels within the UD boundary. 
Stratifications of temperature, LOS velocity, and field strength 
of the UDs were then determined as such averages. The same 
quantities were also determined for the UDs' diffuse background 
(DB), defined as the average over all pixels (ignoring penum- 
bral pixels) in a 400-km-wide ring around the UD boundary. 
UD and DB quantities were retrieved for optical depths between 
log T500 = -2.5 and in steps of 0.5. 

The LOS velocity maps of the inversion result show a clear 
p-mode pattern with a spatial wavelength of about 10" that has 
to be removed to avoid any p-mode influence on our results. The 
usually employed technique of Fourier filtering in 3D few-space 
cannot be applied in our case because only a single map of the 
observational data was available for inversion. We therefore re- 
moved the p-modes in the LOS velocity maps at all used log T500 
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nodes by applying a highpass filter (implemented as the differ- 
ence between the original velocity map and its running boxcar, 
21x21 pixels, filtered counterpart). Since our results depend on 
a careful zero-velocity determination, we re-calibrated the ve- 
locities even if the highpass filter already roughly removed the 
velocity offset. To achieve this we assumed that the dark core of 
the umbra is at rest. The darker part of the umbra is identified 
by thresholding the lowpass-filtered continuum image (11 x 11 
pixels) at 50% of the intensity range. We furthermore excluded 
a circle of 500 km radius around each of the 67 identified UDs 
and subtracted the mean velocity of the remaining dark umbral 
pixels from the velocity maps at each optical depth. This pro- 
cedure was found to be robust in the sense that changing the 
threshold for identifying the darkest part of the umbra by +10%, 
or increasing the radius of the exclusion zone around the UDs by 
200 km, did not influence our results. 

3. Results 

Even at the significantly improved image quality provided by 
the inversion, we were unable to find dark lanes in the cen- 
tral parts of the UDs as rep orted in Schussler & Voglerl (1200 6l) 
in MHD simulations and in iBharti et al.l ( 2007b in other decon- 
volved Hinode images. 

The stratifications of temperature, velocity, and field 
strength, averaged separately over all PUDs and CUDs, are dis- 
played in Fig. |2] While in the upper photosphere (log t^qq = 
-2.5) the considered properties hardly differ between the mean 
UD and DB, they deviate significantly from each other in the 
deep photosphere (log T500 = 0). Compared with their DB, we 
find a temperature enhancement and field- weakening of 610 K 
and 580 G at optical-depth unity for the mean PUD, while for 
CUD the values are 570 K and 530 G. The mean UD magnetic 
field weakens with depth, while the field strength of the mean 
DB increases with depth, as expected. 

The LOS velocity (which is virtually identical to the vertical 
velocity component due to the small heliocentric angle) of DB 
and UD is almost zero in the upper photosphere. Strong upflows 
of -900 m/s and -720 m/s are found at optical-depth unity for 
the mean PUD and CUD, while a weak but significant down- 
flow of 57 m/s is found for the mean DB of the peripheral UDs 
only. These values should be compared with the uncertainty in 
the velocity averaged over the DB of PUDs, cr/ = 10 m/s 
(cr - standard deviation of the 44 DB velocities from their mean 
value, N = 44 - number of UDs). The DB of the central UDs is 
on average at rest within the error bars. A better insight into the 
up- and downflows in and around UDs is given in Figs. [3] and 
|4] The LOS velocity maps at log T500 = are shown in Fig. [3] 
for four typical UDs and reveal that the upflows are mainly con- 
centrated within the marked 200 km vicinity of the UD center, 
while weaker downflows are preferentially found farther out. In 
general, the downflows only partly surround a UD, they are of- 
ten concentrated on one or two sides of the UD. Although weak 
upflows are also present in the vicinities of UDs, the downflows 
dominate for the PUDs, as deduced from Figs. [2] and |4] For a 
better visibility the spatial sampling of the maps was increased 
via bi-linear interpolation. 

In Fig. [4] we averaged the velocities of all pixels within rings 
of 40 km width around the UD center (ignoring penumbral pix- 
els) and plotted the mean velocities as a function of the dis- 
tance to the UD's center, again separately for PUDs and CUDs. 
Note that this method is independent of any determination of 
the UD boundary. Panel (a) of Fig. [4] shows the velocities at 
optical-depth unity. Between - 200 km distance from the UD's 
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Figure 2. Optical-depth dependence of temperature, LOS veloc- 
ity, and magnetic field strength averaged over 44 peripheral um- 
bral dots (top three panels) and 23 central umbral dots (bottom 
three panels). The error bars denote standard deviations of the 
mean (cr/ yN). The red lines exhibit the stratifications of the 
mean UD, while the blue lines correspond to the mean nearby 
diffuse background (see main text for details). Negative LOS ve- 
locity values indicate upflows. 
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Figure 3. LOS velocities at optical-depth unity of typical UDs 
marked as circles of 200 km radius around the position of the 
UD's peak intensity. The UDs in the center of each panel identi- 
fied by letters are the same as in the bottom panel of Fig. Q] 



center, we find upflows (see enlarged velocity maps in Fig. [3] 
where the UDs are marked with circles having a 200 km radius). 
Then, from 200 to 500 km we see downflows (between 200 and 
350 km for CUDs), while for distances larger than 500 km the 
velocity is almost zero. The downflows in the lower photosphere 
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Table 2. Compariso n of UP propertie s at the continuum forma- 
tion height between [Riethmiiller et al.1 (l2008al) and this study. 
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Figure 4. Panel (a)-(c): LOS velocity at constant optical depth 
as function of the distance to the UD center averaged over all 
azimuthal angles of the peripheral (red lines) and central (blue 
lines) UDs. The optical depths are given as text labels. Panel (d): 
Mean continuum intensity profile of the two UD classes. 



peak at a distance of roughly 240 km from the UD's center and 
have values of 1 10 m/s and 58 m/s for the mean PUD and CUD. 
They are thus minute compared to the maximum upfiows in the 
UDs of -960 m/s and -600 m/s. The upfiows and downflows 
are on average stronger for the PUDs than for the CUDs. Both 
up- and downflows increase rapidly with depth (compare pan- 
els (a)-(c)). The upfiows within the UDs are much weaker at 
log T500 = - 1 (panel (b)) and the downflows cannot be seen any- 
more. At log T500 = -2.5 (panel (c)), there is almost no velocity 
signal. The mean intensity profiles are plotted in panel (d) and re- 
veal a half-width-half-maximum radius of 120 km for the PUDs 
and 140 km for the CUDs. 

We next calculated mass fluxes as the sum over all pixels 
within a 500-km vicinity of the UD center since for larger dis- 
tances the velocity is negligible. The required densities were pro- 
vided by the SPINOR code under the assumption of hydrostatic 
equilibrium. Table[TJlists the upward and downward mass fluxes 
per UD, M u , and Md for various optical depths and separately 
for the two UD classes. The uncertainties in Table[T]are the stan- 
dard deviations of the averages over all UDs of a given class 
(cr). The last two columns give the mass flux ratios. The mass 
flux increases strongly with depth due to the density and veloc- 
ity increase. Even if the method of azimuthal averaging leads to 
downflow velocities that are lower than the upflow velocities, the 
averaging is performed over a much larger area so that the up- 
ward and downward mass flux are roughly balanced within the 
uncertainties. 



4. Discussion and conclusions 

We have used a new inversion technique to retrieve the at- 
mospheric parameters of 67 UDs in a sunspot umbra from 
data recorded with the spectropo larimeter onboard Hinode . 
In agreement with earlie r studies (ISocas-Navarro et al ] 12004 
IRiethmiiller et al.Ll2008al) . we found that in the deep photosphere 
the temperature is enhanced and the magnetic field is weakened 
in the UDs compared with their umbral surroundings. Table |2] 
compares the main UD properties retrieved from the conven- 
tional and th e improved invers i on tech nique. For a direct com- 
parison with IRiethmiiller et al.l (I2008al) . who reported peak val- 
ues and not spatial averages, we also listed the peak values ob- 
tained from the new inversion in Table [2] In fact, all values listed 
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in Table [2] are higher for the new inversion method, which em- 
phasizes the considerably improved data quality reached by the 
implicit removal of the telescope's spatial PSF. 

The 2D inversion results revealed clear up flow signals for 
both UD types, while IRiethmiiller et al.l (l2008al) could only find 
them for the PUDs. On average, UDs show upfiows up to a radial 
distance of 200 km from their centers. In general, these upfiows 
are stronger for PUDs than for CUDs. Between 200 and 500 km 
from a UD's center, we found low but significant downflows, 
whereas there is no relevant velocity signal farther away. The 
velocity signal decreases rapidly with atmospheric height. 

Previous observational studies detected upfiows, but 
could not detect downflows associated with UDs (e.g . 



Socas-Navarro et al. , 120041: IRiethmiiller et al 



2008a, 



I Watanabe et al. , 20121) . or at least not systematically ( Ortiz et al. 
l2010l) . This raised the question where all the upflowing plasma 
ends up. The first systematic detection of downflows around 
UDs in this paper gives us the possibility of calculating 
upward and downward mass fluxes. Our finding of the very 
well-balanced mass fluxes depends on the careful velocity 
calibration described in Sect. [2] If all umbral pixels had been 
used for the zero velocity determination, the zero velocity could 
possibly be blueshifted due to the UD upfiows, thus giving rise 
to artificial downflows. This effect is ruled out since our velocity 
re-calibration ignores the UDs and uses the darkest parts of the 
umbra only. 

We furthermore believe that the downflows seen in the top 
left panel of Fig. [4] are real and not a result of ringing effects. 
Such effects could be caused by the nearly axisymmetrical shape 
of the spatial PSF used in our inversion and would have affected 
all quantities. However, plots of temperature and field strength 
versus the UD center distance (not shown) do not exhibi t any 
signs of ringing. According to ISchiissler & Voglerl (120061) . the 
downflows are concentrated in narrow channels preferentially 
at the endpoints of the central dark lanes of the UDs. Our spa- 
tial resolution is insufficient to detect the dark lanes or the nar- 
row downflow channels. The relatively weak downflow signals 
become significantly stronger than the noise only after the az- 
imuthal averaging around UDs. 

The picture introduced in Sc hiissler & Voglerl (120061) of UDs 
as a natural consequence of magnetoconvection in the strong ver- 
tical magnetic field of an umbra is qualitatively confirmed by our 
study. In deep layers the rising hot plasma pushes the field to the 
side, weakening the field there. Around the continuum formation 
height the rising gas cools through radiative losses, turns over, 
and flows down around the UDs. This evidence for overturning 
comes from the mass balance between up- and downflows and 
from the fact that the central upfiows are associated with hot ma- 
terial, whereas the peripheral downflows are cool. Furthermore, 
the rapid decrease of upward mass flux with height is also typi- 
cal for overturning, overshooting convection. The fact that PUDs 
are associated with higher flow velocities and slightly higher 
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Table 1. Upward and downward mass fluxes per umbral dot (computed within 500 km radii) and their ratios at various optical 
depths. 
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1 
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temperature enhancements suggests a more vigorous convection 
than in the central umbra. We note, however, that our results dif- 
fer in one detail from the simulations. Whereas the simulations 
produce concentrated downflows on opposite ends of the UDs, 
the observations reveal a more diffuse downflow structure. 

The 2D inversion technique also returns an enhanced temper- 
ature excess and magnetic -field reduction compared with tradi- 
tional inversions. The temperature and magnetic field anomalies 
turn out to be very similar for PUDs and CUDs (within 10% of 
each other). 

We suggest observations at a spatial resolution higher than 
availa ble here, e.g., with the re-flight of the 1 m Sunrise tele- 
scope dSolanki et al.1 1201 OtlBarfhol et aUl2011l) . or with the en- 
visaged 1.5-m Solar-C telescope, to determine if the UDs show 
the predicted central dark lanes with narrow downflow channels 
at their endpoints. 
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